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Abstract. In a recent article |10| , the authors proved that the non-relativistic Schrodinger 
operator with a generic honeycomb lattice potential has conical (Dirac) points in its dispersion 
surfaces. These conical points occur for quasi-momenta, which are located at the vertices of the 
Brillouin zone, a regular hexagon. In this paper, we study the time-evolution of wave-packets, which 
are spectrally concentrated near such conical points. We prove that the large, but finite, time 
dynamics is governed by the two-dimensional Dirac equations. 
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1. Introduction and Outline. There is great interest within the fundamen- 
tal and applied physics communities in the properties of waves in periodic structures 
with honeycomb lattice symmetry. The (Floquet-Bloch) dispersion relation of such 
structures is known to have conical singularities which occur at the intersections of cer- 
tain bands at high-symmetry quasi-momenta. These conical singularities, also called 
Dirac points or diabolical points, are central to the remarkable electronic properties of 
graphene [151111] and wave-propagation properties in dielectrics (linear and nonlinear) 
with honeycomb structure dielectric parameters [J31 US El [1_7] . Conical points have 
long been known to arise in the dispersion relation of plane waves of the homogeneous 
and anisotropic Maxwell equations [7]. 

In [10] it was proved that for generic honeycomb lattice potentials, V(x), that the 
non-relativistic time-independent Schrodinger equation: 



has conical singularities in its dispersion surfaces. These occur at quasi-momenta 
located at the vertices of the Brillouin zone, 6, a regular hexagon. In this paper we 
prove that the dynamics of solutions of the time-dependent Schrodinger equation: 



for initial data which are spectrally concentrated near the vertices of B, are for very 
large, but finite, times effectively governed by a two-dimensional system of relativistic 
Dirac equations. We next explain our main result, Theorem 15.11 

It is natural to decompose solutions of (|1.2p in terms of its Floquet-Bloch states: 
$ b (x;k)e-^ b ( k )*, where b > 1, k 6 B and n b (k), b > 1 are the eigenvalues of the 
pseudo-periodic eigenvalue problem with quasi-momentum, k; see (|2.9[) - (|2.10[) . At 
a conical singularity (Dirac point) of a honeycomb structure we have two dispersion 
surfaces, graphs of consecutive maps k >-» /x&^k) = jtt-(k) and k >-» /i bl+1 (k) = jtt+(k) 
intersecting conically at each vertex, K* of the Brillouin zone, B: /i* = /i + (K«) = 
/i_(K*). The b\ h and (&i + l) s * spectral bands intersect at the energy fi t and this 
energy is attained by /x_ (k) and fi + (k) at each of the vertices of B. The corresponding 



M = (-A + V(x)) 



(1.1) 



idtt/j = (-A + V(x)) V , 



(1.2) 
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two-dimensional quasi-periodic eigenspace associated with the quasi-momenta K t , 
Nullspace (— A + V — /i*), is spanned by the pair: $!(x;K») and $ 2 ( x ;K«), which 
satisfy the relation: $2(x) = 3>i(— x); see the notion of Dirac point, Definition 13.11 

Theorem [CT] asserts the following for a generic honeycomb lattice potential, V(x): 
Consider initial conditions of the form: 

2 

V(x,0) = 2 «a jO (&0*i(x) , (1.3) 
j'=i 

with fixed, smooth, rapidly decreasing Ojo(X), j = 1,2 and 5 small. We call this a 
wave-packet spectrally localized at K* e B. For such initial conditions the solution 
evolves, approximately, as a slowly modulated superposition of Floquet-Bloch states: 

2 

V(x,t) «e-^ 4 J] ^(Jx^f) $.( x ), (1.4) 

where the modulating amplitudes, Oj(X, T), satisfy the effective Dirac system 

a T a!(X,T) = -Aj (5 Xl +<3 Xa ) aa(X,T) (1.5) 
<5 T o;2(X,T) = -Aj (d Xl -idx 2 ) ai(X,T), (1.6) 

where ^ A« e C. In Theorem 15.11 we establish the validity of (|1.4j) where ol\,o.2 
satisfy (jl.5p - (|1.6l) . on time scales of order (D(S~ 2+e ), for any e > 0. 

To prove Theorem 1 5. 11 we seek a solution of the initial value problem with wave- 
packet initial condition (ll.3[) with leading order term given by the right hand side of 
(|1.4p plus a correction term, rf{x,t), which is represented via the DuHamel formula; 
see (|6.5|) - (|6.7j) . The Dirac equations (|1.5|) - ()1.6[) arise as a non-resonance condition, 
which ensures that rj s (x,t) is small on a time interval: t (D(5~ 2+£ ), for any 
e > 0. Estimation of i] S requires a careful decomposition of the propagator, e -t ( -A+v ')' 
and analysis of its action on functions with quasi-momentum components supported 
near K*, a vertex of B, and those with quasi-momentum components supported away 
from K*. The resonant terms which are removed by imposing equations (| 1 . 5(1 - (jl ,6j) . 
arise from quasi-momenta near K* . A detailed expansion of the normalized Floquet- 
Bloch modes for such quasi-momenta is required. Such modes are discontinuous at 
K* . Components corresponding to quasi-momenta away from K* are controlled, via 
Poisson summation and integration by parts with respect to time, by making use of 
rapid phase oscillations in time. 

Formal derivations of Dirac-type dynamics for honeycomb lattice structures are 
discussed in the physics [15] and applied mathematics [21 [3] literature. A rigorous 
discussion of the tight-binding limit is presented in pQ. Conical singularities have 
long been known to occur in Maxwell equations with constant anisotropic dielectric 
tensor; see, for example, [T3], [7] and references cited therein. 

To put our results in context, we discuss the effective dynamics of two other 
classes of initial conditions: 

1. Ballistic propagation [4 : Take data given by a wave-packet which is local- 
ized at a frequency, p, = /^(K), where ^(k) is regular in a neighborhood 
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of K, p, is a simple eigenvalue of -ff(K) with corresponding Floquet-Bloch 
eigenstate $^(x;K) and VkfJ,r(K) ^ 0: 

Vo(x,0) = 6a (5x) $~ b (x;K) . 
Then, the large time approximate evolution is given by: 
ip(x,t) e"** <5 a(5x,5t) $ 5 (x;K) 

3 r a(X, T) + V k Ms(K) • V x a(X, T) = 0, X = Sx, T = St 

Thus, 

i/>(x, t) * e-** * a ( * • ( x - V k M ? (K) i ) ) $ s (x; K) (1.7) 

for times, i, of order 8~ 2 . 
2. Effective mass (homogenized) Schrodinger evolution [4|: Let K be 

such that /xg(K) occurs at a spectral band (gap) edge. Take wave-packet 

data which is spectrally localized near the frequency /^(K): 

V>o(x,0) = S a (<5x) %(x;K), < <5 « 1. 

Since jt*r(K) is at a band edge, we have Vk^(K) = 0. Furthermore, as- 
sume the Hessian matrix Dj^/i^K) is non-degenerate. Then, the large time 
approximate evolution is given by: 

t/>(x,t) e"^* (5 a(5x,(5 2 t) $ s (x;K) 

where a(X, r) is governed by the constant coefficient Schrodinger equation: 

id r a(X, t) = -V x • A cff V x a(X, r), X = <5x, T = J 2 t (1.8) 

^cff = \ Dln~ b {K) 

for times, t, of the order <5~ 2 . A c g is referred to as the inverse of the effective 
mass tensor. 

1.1. Outline of the paper. In Section[2]we review basic Floquet-Bloch theory 
for general periodic potentials and introduce the class of honeycomb lattice potentials. 
In Section|3]we discuss the main results of the authors' recent paper [10] as well as some 
direct consequences required in the current work. In Section |4] we discuss properties 
of solutions to the two-dimensional Dirac system (|1.5p - (|1.6p . In Section [5] we state 
our main result, Theorem 15.11 on the large, but finite, time Dirac effective dynamics 
for appropriate wave-packet initial data for the time-dependent Schrodinger equation 
with a generic honeycomb lattice potential. The proof of Theorem l5.1l is contained in 
Sections [6] and [7] Appendix [A] gives an elementary proof of the Lipschitz continuity 
of eigenvalues of self-adjoint operators. We thank B. Simon for a sketch of a shorter 
proof using standard perturbation theory; see Chapter XII of |18) . 

In a forthcoming article, we present an analytic perturbation theory of deformed 
honeycomb lattice Hamiltonians, for perturbations which commute with inversion 
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composed with complex conjugation. Conical (Dirac) points persist for small pertur- 
bations of this type, although the conical singularities typically perturb away from the 
vertices of B. These results extend those of [TU] and, in particular, include the case of 
a uniformly strained honeycomb structure. We also consider the analogous question 
of the dynamics of solutions for wave-packet initial data, spectrally concentrated at 
a Dirac point of the deformed honeycomb structure. In this case, the methods of 
the present article apply to establish the large, but finite, time dynamics as being 
given by tilted- Dirac equations. The latter can be mapped to the standard 2D Dirac 
equations by a Galilean change of variables. 

Acknowledgements: CLF was supported by US-NSF Grant DMS-09-01040. MIW 
was supported in part by US-NSF Grant DMS-10-08855. The authors wish to thank 
M. Ablowitz, A.C. Newell and G. Uhlmann for stimulating discussions. 

1.2. Notation. 

1. z e C ==> z denotes the complex conjugate of z. 

2. A, a d x d matrix =^> A 1 is its transpose and A* is its conjugate-transpose. 

3. K m = K mi '" 12 = K + mk = K + miki + m 2 k 2 . 
K,ki and k 2 are defined in Section l2~2l 

4. B denotes the standard Brillouin zone of Figure l2~2l Bh denotes an equivalent 
choice, introduced for convenience in the proofs, centered at K. 

5. V k = e- lkx V x e lkx = V x + ik, A k = V k • V k . 

6. x, y e C", <x, y> = x • y, x • y = xiyi H h x n y n . 

7. For q = q 2 ) e Z 2 , qk = qiki + q 2 k 2 . 

8. <f,g> = ffg 

9. x < y if and only if there exists C > such that x ^ Cy. 

10. We write / = Ox(p) if there exists a constant, C, such that \\f\\x < Cp. 

2. Periodic Potentials and Honeycomb Lattice Potentials. We begin with 
a review of Floquet-Bloch theory of periodic potentials [9) HH [HI [18] . 

2.1. Floquet-Bloch Theory. Let {vi,V2} be a linearly independent set in R 2 . 
Consider the lattice 

A = {mivi + m 2 v 2 : mi, m 2 e Z } = Zvi © Zv 2 . (2-1) 

The fundamental period cell is denoted 

n = { 0i vi + 9 2 v 2 : < 6j sS 1, j = 1, 2 } . (2.2) 

Denote by L 2 er A = L 2 (M 2 /A) , the space of L 2 oc functions which are periodic with 
the respect to the lattice A, or equivalently functions in I? on the torus M 2 /A = T 2 : 

/ e L 2 per A if and only if /(x + v) = /(x), for x e ]R 2 , v e A . 

More generally, we consider functions satisfying a pseudo-periodic boundary condition: 

/ e Ll A if and only if f(x + v) = /(x)e 4k v , for xel 2 , v e A. (2.3) 

We shall suppress the dependence on the period-lattice, A, and write L 2 ., if the choice 
of lattice is clear from context. For / and g in L\ A , fg is locally integrable and A- 
periodic and we define their inner product by: 

</,<?> = f 7Wff(x)dx. (2.4) 
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In a standard way, one can introduce the Sobolev spaces A . 
The dual lattice, A*, is defined to be 

A* = {miki + m 2 k 2 : m 1 ,m 2 e Z} = Zki © Zk 2 , (2.5) 

where ki and k 2 are dual lattice vectors, satisfying the relations: 

ki • Vj = 27T(5,;j . 

If / e L 2 per A then / can be expanded in a Fourier series with Fourier coefficients 

/ = {/m!meZ 2 : 

/(X) = 2 /m e lmk ' X = 2 / mi , m2 e i(-ik 1+ m 2 k 2 ).x _ (2 g) 

meZ 2 (mi,m 2 )eZ 2 

/m ^ pr e-™ k " /(y) dy = ^ e ~^ 1+m ^ 2 y y /(y) dy (2 J) 

Let V^x) denote a real- valued potential which is periodic with respect to A, i.e. 

l/(x+v) = y(x), forxeM 2 , veA. 

Throughout this paper we shall also assume that 

V e C X '(R 2 /A) . (2.8) 

We expect that this smoothness assumption can be relaxed considerably without much 
extra work. 

For each k e R 2 we consider the Floquet-Bloch eigenvalue problem 

H v $(x;k) = /x(k) $(x;k), xeM 2 , (2.9) 

$(x + v;k) = e jk v $(x;k), v e A, (2.10) 

where 

H v = -A + V(x) . (2.11) 

An L\r solution of (I2.9[) - (|2.10[) is called a Floquet-Bloch state. A function which 
satisfies the boundary condition (|2 . 10[) is said to be k— pseudo-periodic. 

Since the eigenvalue problems (|2.9p - (l2.10l) are invariant under the change k i— ► 
k + k, where k e A*, the dual period lattice, the eigenvalues and eigenfunctions of 
(|2.9[) - (|2.10j) can be regarded as A*— periodic functions of k, or functions on R 2 /A*. 
Therefore, it suffices to restrict our attention to k varying over any primitive cell. It 
is standard to work with the first Brillouin zone, B, the closure of the set of points 
k G M 2 , which are closer to the origin than to any other lattice point. 

An alternative formulation is obtained as follows. For every k e B we express the 
Floquet-Bloch mode, $(a;;k), in the form 

$(x; k) = e lk x p(x; k) . (2.12) 

Then p(x; k) satisfies the periodic elliptic boundary value problem: 



H v (k)p(x;k) = /x(k)p(x;k), xeR 2 , 
p(x + v;k) = p(x;k), v e A, 



(2.13) 
(2.14) 
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where 

H v {k) = - (V + ikf + V(x) . (2.15) 

The eigenvalue problem (|2 .9[) - (|2 . 10[) . or equivalently (|2.13[) - (|2.14[) . has a discrete spec- 
trum: 

M i(k) < M2 (k) < M3 (k) < ... (2.16) 

with eigenpairs p&(x; k), /ib(k) : 6 = 1, 2, 3, ... . The set {pb(x; k)}fc^i can be taken 
to be a complete orthonormal set in L 2 er (R 2 /A). 

The functions /i&(k) are called band dispersion functions. Some general results 
on their regularity appear in [191 [5] . 

Since V is assumed to be smooth, elliptic regularity theory implies for each b > 1 
and ke8, that 

x^ Pb ( X ;k) is C"°(E 2 /A). 
Furthermore, there exists a constant C& ^y, depending only on 6,/3 and V", such that 

max ||(5^ 6 (-;k)|| LK , (n) < C b ^y. (2.17) 



We shall also require the regularity of the mapping k i— > /if,(k). 

Proposition 2.1. The eigenvalue maps k >— > /Xb(k), b ^ 1, are Lipschitz contin- 
uous. 

Proposition ^. II (see also Proposition lA.2j) is a consequence of the general result on 
Lipschitz continuity of eigenvalues of second order elliptic operators (Theorem lA.lj) . 
stated and proved in Appendix [A] 

Remark 2.1. Although the eigenvalue maps, k >— > /ib(k) , are Lipschitz functions 
the Floquet-Bloch mode maps, k i— > Pb(x; k), are in general not even continuous \1 9f . 
Indeed, we shall see this behavior explicitly in a neighborhood of degenerate eigenval- 
ues; see Theorem \3. 2\ 

As k varies over B, /Xb(k) sweeps out a closed real interval. The spectrum of 
—A + V(x) in L 2 (R 2 ) is the union of these closed intervals: 

spec(iiy) = (J spec(ffy(k)) . (2.18) 

Moreover, the set 

|J |J {$ b (x; k)}, $ 6 (x; k) ee e jk x p b (x; k), (2.19) 

b»l kei3 

suitably normalized, is complete in L 2 (M. 2 ): 

/(x) = J f <$ b (.;k),/(-)> L2(R2) $b(x;k) dk = £ / A« ^k) dk , 

(2.20) 

where the sum converges in the L 2 norm. 
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Moreover we have, with respect to the Floquet-Bloch basis, the Plancherel The- 
orem: 

\\f\\h m = £ / i^' 2 dk ( 2 - 21 ) 



Remark 2.2. The A*— •periodicity of the Floquet-Bloch modes implies that we 
can express (|2.20p equivalently in terms of a dk— integral over any fundamental period 
cell. A convenient choice, to be used below, is one where the integal over B is replaced 
by an integral over 

B h = K + B . (2.22) 

That K is an interior point to this fundamental domain, rather than a vertex, will 
simplify certain computations below. 

Thus it is natural to introduce Sobelev spaces, defined in terms of the Floquet- 
Bloch coefficients as follows: 

II/IW 2 ) * IIU+|tf| 2 ) f /lli 2 ( R2 ) = 2 I (1 + Mk)| 2 )* |/ b (k)| 2 dk 

*£(1 + |&| 2 )« / l/Uk)| 2 dk. (2.23) 

63*1 J B 

The latter approximation is a consequence of: 

| W (k)|~|6|, (2.24) 

The Weyl law (I2.24p holds uniformly in B. 

Note the simple consequence of (|2.23|) . to be used in Section [TJ 

fb = 0, for b outside a fixed finite set => ||/||j?»(r2) < ||/||z, 2 (k 2 )! s ^ 0. (2.25) 

2.2. The period lattice, A h , and its dual, A£. 

Consider A/, = Zvi ©Zv 2 , the lattice generated by the basis vectors: 



vi = a " , v 2 = a " , a > 0. (2.26) 















(!) 


| , v 2 = a | 

















Note: Ah ("h" for honeycomb) is a triangular lattice, that arises naturally in connec- 
tion with honeycomb structures; see Figure [2TT1 

The dual lattice A* = Zki © Zk 2 is spanned by the dual basis vectors: 



1 \ / l 

2 \ / 2 



4tt 



2 / \ 2 
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1.5- 



1- 



0.5- 



-0.5- 





-1.5- 



-1.5 



-0.5 0.5 



1.5 



Fig. 2.1. Part of the honeycomb structure, H. H is the union of two sub-lattices Aa = A + A^ 
(blue) and Ab = B + (green). The lattice vectors {vi, V2} generate A^. 



where 



k e ■ ve = 2nSw , 

|vi| = |v 2 | = a, vi • v 2 

|ki| = |k 2 l = g, ki-ka 



a 
~2 



The Brillouin zone, B, is a hexagon in M 2 ; see Figure 
the vertices of B given by: 

K=i(ki-k a ), K' = -K = i (k a - k x ) 



(2.28) 
(2.29) 
(2.30) 

Denote by K and K' 
(2.31) 



All six vertices of B can be generated by application of the rotation matrix, R, which 
rotates a vector in R 2 clockwise by 27r/3. R is given by 



R 



2 

VI 

2 



v3 

2 

_ 1 

2 



(2.32) 



and the vertices of B fall into two groups, generated by action of R on K and K': 

K type - points : K, RK = K + k 2 , R 2 K = K - k x 

K' type - points : K', RK' = K' - k 2 , R 2 K' = K' + ki . (2.33) 



Charles L. Fefferman and Michael I. Weinstein 



9 



2r 

1.5- 



1- 




-1- 



-1.5- 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 



Fig. 2.2. Brillouin zone, B, and dual basis {ki,k2}. K and K' are labeled. Other vertices of 
Bh obtained via application of R, rotation by 2ir/3; see equation l |2.33|l . 



Remark 2.3 (Symmetry Reduction). Let ( $(x;k),^(k) ) denote a Floquet- 
Bloch eigenpair for the eigenvalue problem (|2.9p - (|2.10p with quasi-momentum k. Since 
V is real, ( <I>(x;k) = $(x;k),/i(k) ) is a Floquet-Bloch eigenpair for the eigenvalue 
problem with quasi-momentum — k. Recall the relations (|2.33[) and the A^- periodicity 
of: k i— * /i(k) and k i— > <£>(x;k). It follows that the local character of the dispersion 
surfaces in a neighborhood of any vertex of B is determined by its character about any 
other vertex of B. 

We reminder the reader that below, it will be convenient to work with Bh = K + B 
as our Brillouin zone, explained in Remark 12.21 

2.3. Honeycomb lattice potentials. 

For any function /, defined on R 2 , introduce 

ft[/](x) = /(i?*x), (2.34) 

where R is the 2x2 rotation matrix displayed in (|2.32[) . 
Definition 2.2 (Honeycomb lattice potentials). 

Let V be real-valued and V 6 C X (IR 2 ). V is a honeycomb lattice potential if there 
exists Xo e R 2 such that V(jt) = V^(x — Xo) has the following properties: 

1. V is Ah— periodic, i.e. V(x + v) = V(k) for all x e R 2 and ve Aj. 
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2. V is even or inversion-symmetric, i.e. V(— x) = V(x). 

3. V is TZ- invariant, i.e. 

TZ[V](x) = V^(i?*x) = V(x), 

where, R* is the counter-clockwise rotation matrix by 2ir/3, i.e. R* = R , 
where R is given by (|2.32l) . 

Thus, a honeycomb lattice potential is smooth, A^- periodic and, with respect to 
some origin of coordinates, both inversion symmetric and TZ- invariant. 

Remark 2.4. As the spectral properties are independent of translation of the 
potential we shall assume in the proofs, without any loss of generality, that xq = 0. 

Remark 2.5. A consequence of a honeycomb lattice potential being real-valued 
and even is that if ($(x; k), fi) is an eigenpair with quasimomentum k of the Floquet- 

Bloch eigenvalue problem, then {^{— x; k), fi^j is also an eigenpair with quasimomen- 
tum k. 

A key property of honeycomb lattice potentials, Vh, used in our spectral analysis 
of — A + Vh [TO], is that if K* denotes any vertex of Bh, then we have the commutation 
relation: 

[H, Hv h (Ki,)] = . (2.35) 

It is therefore natural to the split L 2 ^ , the space of K»- pseudo-periodic functions, 
into the direct sum: 

where L 2 ^ g are invariant eigen-subspaces of TZ, i.e. for a = l,r, r, where r = 
exp(27ri/3), and 

4„, = {geL 2 K ^.TZg = ag} . (2.37) 

3. Dirac points. 

We begin with a precise definition of a Dirac point. 

Definition 3.1. Let V(x) be a smooth, real-valued, even (inversion symmetric) 
and periodic potential on M 2 . Denote by B the Brillouin zone given in Remark \2.2\ 
We call K e B a Dirac point if the following holds: There exist an integer b\ ^ 1, a 
real number and strictly positive numbers, A and S, such that: 

1. /i* is a degenerate eigenvalue of H with K pseudo-periodic boundary condi- 
tions. 

2. dim Nullspace(^H — = 2 

3. Nullspace^H — fj^l) = span|$i(x), $ 2 (x)|, where <I>i e L 2 ^ T and $ 2 (x) = 

4- There exist Lipschitz functions /i+(k), 

Mbi( k ) = M-( k ) Mbi + i( k ) = M+( k ), M±( K ) = M* 

and E+(k), defined for |k — K| < S, and k— pseudo-periodic eigenf unctions 
of H: $+(x;k), with corresponding eigenvalues /i+(k) such that 

/i+(k)-/i* = + A |k — K| (1 + E+Ql) ) and 

M-(k)-M* = - A |k-K| ( 1 + E-Ql) ), (3.1) 

where \E+(k)\ < C|k-K| for some C > 0. 
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Remark 3.1. In UOj/ we prove the following 

Proposition 3.2. Suppose conditions l.,2. and 3. of Definition COl hold and 
denote by {c(m)} m( =s the sequence of T Fourier- coefficients o/<I>i(x). Define the 
sum 

a„ = 2 c ( m ) 2 ( !V K ™ ' (3 - 2) 

with S cz Z 2 ; see \101. If A| ^ 0, i/ien 4. of Definition HOI /to/ds ('see dSHJ) j ure'i/i 
A=|A„|. 



We next recall the statement of Theorem 5.1 of pj)| concerning the existence of 
Dirac points for the Schrodinger operator with a generic honeycomb lattice potentials. 

Theorem 3.1. Let Vh(x) honeycomb lattice potential. Assume further that the 
Fourier coefficient ofVh, V\,i, is non-vanishing, i.e. 

V hl = f e-^+^-y V h (y) dy # . (3.3) 
Jn 

Consider the one-parameter family of honeycomb Schrodinger operators defined by: 

ff« = -A + e y h (x) . (3.4) 

There exists a countable and closed set C cz R suc/i that for all e ^ C, the vertices, K*, 
o/Sh are Dirac points in the sense of Definition \3.1\ 

More specifically, the following holds for e ^ C: There exists b\ ^ 1 such that 
/i* = /i^ (K*) = /ij +1 (K*) is a K* pseudo-periodic eigenvalue of multiplicity two 
where 

1. /il is anL\ Lr - eigenvalue of of multiplicity one, and corresponding 
eigenf unction, $f(x). 

/il is an f - eigenvalue of of multiplicity one, with corresponding 

eigenf unction, $|(x) = x). 

\i% is not an x - eigenvalue of H^K 

2. There exist S e > 0, C c > and Floquet-Bloch eigenpairs: ($^_(x; k), /i^(k)) 
and ($l(x; k), /il(k)) 7 and Lipschitz continuous functions, E+(h), defined 
for |k — K*| < S e , such that 

/4(k)-//(K*) = + |A|| |k-K,| (1 + E%(k) ) and 
M ;(k)-//(K,) = -|A|| |k-K,| (l + Sl(k)), 

w/iere 

A§ = X! c ( m ^ £ ' e ) 2 ( n K ? * ( 3 - 5 ) 

is given in terms o/{c(m, /x e , e)}, £/ie i 2 ^ T - Fourier coefficients o/$ e (x; K*). 
Furthermore, |i?+(k)| ^ C e |k — K*|. Thus, in a neighborhood of the point 
(k, n) = (K*,/l«J) £ M 3 , the dispersion surface is closely approximated by a 
circular cone. 
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3. There exists e° > 0, such that for all e e (— e°, e°)\{0} 

(i) eVi i > ==> conical intersection of I s * and 2 nd dispersion surfaces 

(ii) eVi t i < ==> conical intersection of 2 nd and 3 rd dispersion surfaces . 

Our point of departure in this paper will be a periodic Schrodinger operator, 
—A + V(x), where V is a honeycomb lattice potential. We fix a Dirac point, ensured 
to exist by Theorem 13. 11 and study the large time dynamics of wave-packets which 
are initially (for t = 0) spectrally localized near these points. Thus we have two band 
dispersion surfaces, ksU^ M&iO*) = M-O 4 ) an d k >-* Mbi+i(k) = M+(k) which touch 
conically at k = K» with /i+(K*) = 

Let 

$i(x) = dc f $+(x) e r 

$ 2 (x) = $i(-x) = dcf $_(x) e f (3.6) 

span the two-dimensional subspace of the degenerate eigenvalue [i* for H: 

H*j(x) = (-A + V)* s = frfy (3.7) 
$ 3 -(x + v) = e lK * v $j-(x), xeR 2 , vsAj (3.8) 

We also define the periodic vectors: 

pi(x) = e-***-*^*), p 2 (x) = e- iK * x $ 2 (x) . (3.9) 

We choose these states to be orthonormal 

<$z,5> m > L2(n ) = Si m , i,m=l,2. 



To study the time evolution (|1.2j) we expand the solution of the initial value 
problem with data ipo using the complete set of Floquet-Bloch modes: 

^ lHt ^>o = Yi[ e-^( k ) t <$ b (-;k),^ (-)><i> b (x;k)dk (3.10) 

Suppose H has a Dirac point, K. We remark that in ()3.10p we choose the Brillouin 
zone, Bh = B+K., which is centered at K; see Remark [2. 2 1 For initial conditions which 
are spectrally supported near K, the time evolution depends on the precise behavior 
of the Floquet-Bloch modes $+(x;k) for k near K. While Theorem 13.11 shows that 
the eigenvalues ^+(k) are Lipschitz functions in a small neighborhood of k = K, the 
eigenfunctions $+(x; k) are not continuous functions of k in a neighborhood of K. 

Theorem 3.2. Let K denote a Dirac point in the sense of Definition \3.1\ In 
particular, let ($+(x; k), /i+(k)) denote the k— pseudo-periodic eigenpairs as in part 
4 of Definition \3. R Introduce the A^— periodic functions p+(x;k) by 

$+(x;k) = e 4k x p+(x;k), <p (. ; k),p 6 (.; k)> = S ab , a, b e {+, -} . (3.11) 

Let k = K + k. Then, for < |«| < S and pj given by (|3.9p we have: 



(i+(K+k) = fi* ± | Ao | (k 2 + k 2 )2 + C>(k 2 + k 2 ), 



(3.12) 
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and p+(x; k), a priori defined up to an arbitrary (complex) multiplicative constant of 
absolute value 1, can be chosen so that: 

1 K\ + i&2 1 
p ± (x;K + .) = -_^_^( X ) ± _p 2 (x) 



+ 0^( R2 /A h) ((^ + «2) i ) • (3-13) 



Proof of Theorem \3.2[ The proof builds on the proof of Theorem 4.1 of [lOj . in which 
the k— pseudo-periodic eigenvalues, /Li+ (k) , are constructed. These were shown to 
be Lipschitz continuous functions for k varying in a neighborhood of K. We now 
consider the associated Floquet-Bloch modes. 

Let k— pseudo-periodic Floquet-Bloch modes can be expressed in the form 

$(x;k) =e lkx p(x;k), 

where p(x; k) is Ah periodic. Since we are interested in the character of p(x; k) for 
k near K we set: k = K + k. Then, p(x; K + k) satisfies the periodic eigenvalue 
problem: 

H(K + k) p(x; K + k) = fi(K + k) p(x; K + k) , (3.14) 
p(x + v; K + k) = p(x; K + k), for all v e A h , (3.15) 

where 

HQs) = - (V x + «k) 2 + V(x) . 

The eigenvalue problem p.l4j) - (|3.15|) has eigenvalues, computed via degenerate per- 
turbation theory of the double eigenvalue /i* of iJ(K), given by: 

/i+(K + k) = ^ + /i« (3.16) 
= ±|A„| | K | + 0{\n\ 2 ) . (3.17) 

Denote by Q± the projection onto the orthogonal complement of span{pi,p2}- 
Then, 

i? K (M = (H(K) -^iy 1 : Q ± L 2 (R 2 /A h )^Q ± L 2 (R 2 /A h ) 

is bounded. Furthermore, via Lyapunov-Schmidt reduction of the eigenvalue problem 
f|3. 14|) - (|3. 15|) we obtain, corresponding to the eigenvalues /x+(K + k) the Floquet-Bloch 
modes: 

p+(x;K + K ) = (/ + R K (fu)Q± (2t/s • (V + *K)) ) ( ap x (x) + /3p 2 (x) ) 

+ O ff2(R2/Afe) (| K |(H 2 + |/3| 2 )^) . (3.18) 
The k— dependent coefficients a and /3 satisfy the homogeneous equation 

7W(/i (1) ,k) ( 1 = 0, (3.19) 
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where M.{^\ n) is a 2 x 2 matrix of the form: 

/ M W +C(|k| 2 ) -^(«i +m 2 ) +0(|k| 2 ) 

V -Aj(«i - feu) + (|k| 2 ) ?W + O (N 2 ) 

see [10] for the derivation of ([3l8 ]) -([3l9 ]l , 

The Floquet-Bloch modes, p+(x;k), are finally determined by the solutions a = 

a+(/i), j3 = P+{k)) of the homogeneous algebraic system f|3 . 19[) for the choices /z+'(/c) 
in p. 171) . We select (normalized) solutions of (|3.10[) as follows: 



^ ] = ±\\\ H + C%| 2 ), 



1 K1+iK +0 / (k 2 + k 2)1 

±71 + O ((«? + *§)* 



(3.20) 



Here, we take advantage of the observation that for a rank-1 matrix 

B 



A 



lies in the nullspace of A. 



A B 
C D 



the vector 



-A 



Thus, with k = k — K, we have upon substitution of (|3.20[) into p,18j) the 
expansions of fi+ (k) and p+ (x; k) in p,12|) and (|3. 13|) . 



4. 2D Dirac equation. In this section we collect results on well-posedness and 
estimates on solutions of the two-dimensional Dirac system (|1.5I) - (I1.6[) . Taking the 
Fourier transform of (|1.5[) - (|1.6[l wc obtain for a(H, T) = a(£i,£2,T) the equation 



idn 



0(H) 



0(H) 



Oil 

a 2 



, where 



M& + i£ 2 ) 



0(H)* 



(4.1) 



Remark 4.1. It follows from (|4.ip £/ia£ £/ie dispersion relation for (|1.5I) - (|1.6[) is: 



f/ie effective dynamics are non-dispersive. 

Remark 4.2. The system (|1.5[) - (|1.6[) /ias £/ie structure of Dirac system: 



idi 



Oil 
Q 2 



1 d 



1 <9 



i 3Xi 



dX 2 



Oil 

a 2 



(4.2) 



T7ie 2x2 matrices a\ and o~\ satisfy the relations 



lAjj 2 Id, and a\a\ + o\o\ 



Here, Oj, j = 1,2 are given by the Hermitean matrices: 

a\ = Ajcti, <r| = A b ct 2 



(4.3) 
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in terms of Aj and the standard Pauli matrices o~\ and er 2 •' 

*)•--(! «) <«) 

Finally, note that ctj, j = 1,2 satisfy the two-dimensional wave equation, with wave- 
speed |A(j|: 

=|A„| + 

Proposition 4.1. 
ylsswme a(X,0) = a (X) e (# S (R 2 )) 2 . T/ien, 

d%a e i rx ([0,oo) ; (i7 s - fc (M 2 )) 2 ) , /or k < s . 

In particular, 

1. The Fourier transform of the solution, a(X, T) is given explicitly by 

a(S,T) = e-* n ( H > T d (H), 0(H)* = fi(H) 

2. ForaHBeM 2 , |d(S,T)| = |a (H)| and therefore 

3. For any aeZ 2 with |a| < s 

|| <%a(X,T) || i2(R2) = || 5fta(X,0) || L2(R2) (4.5) 



5. Effective Dirac dynamics; statement of main result, Theorem l5.ll A 

general solution of the time-dependent Schrodinger equation, constrained to the de- 
generate 2-dimensional eigenspace associated with eigenvalue, /z* = /i(K*), associated 
with the Dirac point, K, , is of the form 

V>(x,t) = e"^* ( ax $i(x) + a 2 $ 2 (x) ) , (5.1) 

where ct\ and a 2 are arbitrary constants. 

Consider now a wave packet initial condition, which is spectrally concentrated 
near K*: 

^(x,0) = VoW = S (a 10 (<5x) $ 1 (x)+a 20 ((5x) $ 2 (x) ) 

= S ( a 10 (5x) pi(x) + a 20 (<5x) p 2 (x) ) e iK *' x (5.2) 

Here, (5 is a small parameter. We assume aio(X) and a 2 o(X) are Schwartz functions 
of X. We expect that this assumption can be weakend considerably without difficulty. 
The overall factor of 5 in (|5.2[) is not essential (the problem is linear), but is inserted 
so that tpQ has L 2 (M. 2 )- norm of order of magnitude one. 
We seek solutions of (|1.2[) . (15.21) in the form: 

ip s (x,t) = e - *^' ^ 5 oy(<5x,<ft)$,,(x) +?7 5 (x,i)j . (5.3) 
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where ^(x, 0) = 0, a 3 (X, 0) = a j0 (X), j = 1, 2 (5.4) 
to ensure the initial condition (|5.2j) . 

The goal is to show that the Schrodinger equation (|1.2p has a solution of the form 
(|5.3p with an error term, n s (x,t), which satisfies 

sup || V \-,t) ||^(k 2 ) = 0(5 T *), 5^0. (5.5) 

for some r* > 0, provided the slowly varying amplitudes Oj(<5x, St), j = 1,2 evolve 
according to the system of Dirac- type equations (|1.5I) - (|1.6[) . Here, p > and E\ > 
are arbitrary. 

We shall prove the following 

Theorem 5.1. Assume 



a (X) = 




and let d?(X,T) denote the global-in-time solution of the Dirac system (|1.5[) - (|1.6p with 
initial data d?(X,0) = d?o(X). Consider the time- dependent Schrodinger equation, 
(|1.2[) . where V(x) denotes a potential for which the conclusions of Theorem \3.1\ hold, 
e.g. V(~x) = eVh(x), where Vh is a honeycomb lattice potential satisfying V\_\ ^ 
and e is not in the countable closed set C. Assume initial conditions, tpo, of the form 
(|5.2p . Fix any p > and S\ > 0. Also choose non-negative N = {n\,n2) G Z 2 with 
n\,ni 5= 0. Then, (|1.2p has a unique solution of the form (|5.3p . where for any \a\ < N 

sup ||3^(x,t)|| = o(<T*), 5^0 (5.6) 

for some t* > 0. 

6. Proof of Theorem 15.11 We begin with a summary of consequences of The- 
orems 13.11 and 13.21 which we use in the proof. 

Recall that the spectral bands p + (Bh) and p-(Bh) touch conically at the the ver- 
tices of Bh- Specifically, p+(k) = M-(k) = p* for k = K* equal to any of the three K- 
type vertices, {K, RK, R 2 K} and any of the three K'- type vertices, {K', RK' , R 2 K'}, 
where K' = — K. Here, R denotes the 7r/3- counterclockwise rotation. 

For a Dirac point K*, equal to any vertex of Bh' 

(PI) There exist L 2 ,- eigenvalues of H, denoted p+(k), such that 

for all k e Bh satisfying |k — K*| < kq 

p+(k) - M (K,) = + |A„| |k-K,| (1 + E+(k — K*)) , (6.1) 

where E+(n) are Lipschitz continuous in k and E+(n) ^ C+ \k\. Here, Aj is 
a constant given by (|3.5p . 
(P2) There are constants k,\,C\ such that 



for all k £ Bh satisfying |k — K,| < n\ and all b $ { + , — } 
|/i 6 (k)-/i(K*)|>Ci>0. (6.2) 
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(P3) For < \k\ < 5 sufficiently small 

P+(x;K* + k) = -^(a(K)pi(x) + p 2 (x) ) + O ff 2 (R 2 /A(i) (|k|) , (6.3) 

where «(«;) = (6.4) 

(P2) is a consequence of (PI) and the continuity of the eigenvalues k >-» /i;,(k;); see 
Proposition ^. II 



Without loss of generality and for simplicity: 

1. we take 

K, = K 

2. we recall that the Brillouin zone, Bh, is centered at K; see Remark 1 2. 2 1 Since 
we assume wave-packet initial data which are spectrally localized near K in 
M 2 /A^, this equivalent choice, which puts K on the interior of the Brillouin 
zone, simplifies the analysis. 

To prove Theorem 15. 11 we study the evolution equation for 77 (x, t) obtained by sub- 
stitution of (1531) into (fl~2"j) : 

id t r] s = (H-fM^V 5 
2 



S 2 



J] [ id Taj (X,T) $,(x) + 2V X ^-(X,T) • V x $,(x) 



2 



(X,T) = ((5x,(5t) 



„6 



5 3 V A x a 3 -(X,T)$ 3 -(x) , (6.5) 

J (X,T)=(5x,«) 



^(x,t = 0) = . (6.6) 

6.1. Estimation of the error, -q 6 (x,t). Using the DuHamel principle we may 
rewrite (|6.5[) as the equivalent integral equation: 

jj 4 (.,t) = i5 2 V f e -^-^)(t-s) [ idsa] (S;8s)$ 3 (-) + 2Vxa ] (5;5s)-X7 x $ J (-)]ds 

2 

+ id 3 e-* (fl, -"* )Ct — } Axcy(M«)*i(-) ds (6.7) 

3=1 ■'0 

The second integral in (I6.7P (call it /nt2) can be bounded as follows. Let s denote 
an even positive integer. Recall from elliptic theory that 

WfWh < \\f\\ 2 + |(-A + y + c)*/l3 s||/llr. , (6.8) 

where c > sup xeR 2 |V(x)|. Using (|6.8j) . we can bound ||/nt2||ff s in terms of ||/n<2||L 2 
and || (—A + V + c)^ Int^ih 2 ■ Taking the L 2 norm of Int 2 and using that e~ tHt is 
unitary in L 2 , we obtain 



2 nt 

Int 2 \\ L 2 <S 3 Y lAx. aj {S',S8)lv ds - 
3=1 Jo 



(6.9) 
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Next, using that H = —A + V commutes with e~ lHt we obtain 

2 r t 



\\(-A+V + c)ilnt 2 \\ L 2 <5 3 Y / \\(-A+V + c)iA xaj (S-,Ss)\\ L 2 ds. (6.10) 

3 = 1 J ° 



Next note, via Proposition l4.il that 

ll a j(^-!^)ll// m (R 2 ) < |ao||ff™(R2) 
Let ei > be arbitrary. Then, equations (|6 . 1 1[) imply 



(6.11) 



sup || in£2 1| ir s 



sup 



2 

<5 3 Y / e i(K- M *)(*-«) A x ai(*-,<5s) $ jO ds 

.,=1 ./o 



(6.12) 



It therefore suffices to estimate the first time-integral in (|6.7[) . This time- integral 
is the solution of the initial value problem of the form: 

id t f{*,t) - (H-^)f s (x,t) = 5 2 £ 7r («5 X ,<ft)* r (x), (6.13) 

r 

/ 5 (x,0) = 0, (6.14) 
where we find it convenient to introduce the notation 

2 

Vj 7r (X,T)* r (x) = - J] [ id T aj(X,T) $i(x) + 2V x ai(X,T).V x $ J (x) 



(6.15) 



where 



7i(X,T) = i^a^X.T), * 1 (x) = $ 1 (x), 
72 (X,T) = id T a 2 (X,T), * 2 (x) = $ 2 (x), 



(6.16) 
(6.17) 

73 (X,T) = 2V x «i(X,T), * 3 (x) = (*3 I i(x),*3,2(x)) = V x 4i(x) 1 (6.18) 
74 (X,T) = 2V x a 2 (X,T), * 4 (x) = (^i(x),f 4 , 2 (x)) = V x $ 2 (x) . (6.19) 



Note that *j(x) e L^, 1 j ^ 4. 

By hypotheses on a.j and $j, and Proposition l4.il the functions 7r and ^> r satisfy 
the following properties: 



(6.20) 
(6.21) 



7r (X,T) is C°°(M T ;5(K X )) 

||7r(-,r)| H . (Ra) < ||a(-,r = o)| H . + i 

¥ r (x) e C°°(M X ) n 

and therefore satisfies the pseudo-periodic boundary condition 

* r ( x + v) = e iK v * r (x) . (6.22) 



We also write ^f r in the following useful form: 

¥ r (x) = e jKx 7> r (x), P,.(x + v) = V r (y) , v e A h 



(6.23) 
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By (glTl), for any/3eN2, 

Wd^Vrh^) < C r ,/3 • (6.24) 
Theorem 15. II is therefore reduced to the following 

Proposition 6.1. Let ai(X,T) cmdo^X, T) satisfy the system of Dirac equa- 
tions: 

d T a x = — Ay (d Xl + idx 2 ) «2 (6.25) 

d T a 2 = —Xi(dx 1 —idx 2 )ai (6.26) 

with initial conditions as in Theorem 15.11 Then, for any s ^ l,p > and Ex > 0, 
there exists a unique solution, f s (x,t) of (|6.13|) . which satisfies the estimate 

sup || f\-,t) || < C <5^, for (5 1 . (6.27) 



So to finish the proof Theorem 1 5. 1[ we only require a proof of Proposition ^. II 
7. Proof of Proposition 16.11 By the completeness of Bloch modes, 



V / e -*(w(k)-M.)t <$6 (. ; k), 5 (.)> L2(R 2) $ b (x;k) dk (7.1) 



Thus, 



f(x,t)=Y J [ / b 5 (k,t)$ b (x;k) dk, (7.2) 

where 

#(M) = -i6 2 f dse-^M-^-*) /$ 6 (.;k),2 7 r(*-,*«)*r(-) ) • (7-3) 

\ r / L 2 (R 2 ) 

VFe s/ia/Z henceforth omit the superscript S from f 5 and f£ . 

Decompose / into /u, frequency components which lie in the two spectral bands: 
fi + (Bh) and fi-(Bh) intersecting at the Dirac point /i*, and fn c , frequency components 
which lie in all other spectral bands: 

/(x,t) = /o(x,t) + f D o(x,t), where 

fn(x,t) = J] / / 6 (k,t)$ 6 (x;k) dk (7.4) 



be{ + ,-} 



and 



f D c( X ,t) = 2 / / b (k,t)$ b (x;k) dk. (7.5) 

b<f{ + --} JBh 

Let's focus initially on fu(jc,t). We distinguish between frequencies which are 
"near" and "bounded away from" /i*, as these correspond to whether the complex 
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phase in (|7.1[) is non-oscillatory or, respectively, oscillatory in t. Recalling property 
(P2) and our choice of Brillouin zone with K in its interior: we further decompose fu 
into its quasi- momentum components near and away from from any of the points K: 

/ B (x,t) = J] f X (|k-K| <6 r ) / 6 (k,i) $ b (x;k) dk 

be{ + ,-} J ®» V ' 

+ S / x(^ r < |k-K|) / b (k,t) $ fc (x;k) dk 

be{ + .-} jB h 

S3 ^ / [//^(k,t) +/j/, 6 (k,t)] $ 6 (x;k)dk 
= 2 ( /j, 4 (x,i) + ///, 6 (x,i) ) 

= // >D (x,i) + f n , D (x,t) ■ ( 7 - 6 ) 
Here, < r < 1 will chosen less than but close to 1. By (|2.21[) and (|2.23l) 

!/('>*) lif«(R2) * I/d('> *)lli=(M 2 ) + ll/-D c ('' *)liT»(R2) 

* 2 IIA-d(''*)!Il 2 (r 2 ) + ll/D c (''<)llr»(E 2 ) 
,/=/,// 



(7.7) 



J=IJI b=-_ 



where we have used (|2.25p . We show below, for fixed p > 0, e± > that each term in 
(1771) is 0(S £1 ) for < t < p (T 2 + £l as 5 I 0. 

In the calculations below we shall require a detailed expansion of inner products 
of the form: 



<$ 6 (. ; k),r(*.,*«)*(.)> £aCHa) , 

where T = r(X, T) is in Schwartz class and e L^, i.e. 

*(x) = e iK - x P(x), 7>(x + v)=7>(x), veA h ; 
see (|7.3p . The following proposition will be used: 



(7.8) 



(7.9) 



Proposition 7.1. Let T(X,T) denote a Schwartz class function o/X, varying 
smoothly in T . Denote by T(5, T) i£s Fourier transform with respect to the X variable. 
Then, 



<$ b (-;k),r(<Ms)*(-)W) 

-2 y e *m- y p i k, + //) jk. + i k - K i 



/ My; k) 

Jn 



Ss 



V{y) dy 
(7.10) 

Proof of Proposition^^ Recall from ([2~TO1) and (joT2"3")) that $ b (x; k) = e ikx p 6 (x; k), 
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where Pf,(x + v; k) = Pfc(x; k) for any v e A/,. Thus, using (|7.9p . we find that 

<* b (.;k),r(5-,^)*(.)> i2(R2) 

= / e-^ K )yr(^y,fc)^(yTkTP(y)dy 

= V / e -«(k-K).(y+mv) p^ y + ^ pb ( y + mv; k ) V (y + mv ) dy 

mez 2 •'n 

= S / e- i(k - KHy+mv) r(J(y + mv),J S )^(yTkT7'(y) dy 

meZ 2 , ' n 

= / | 2 e- j(k - KHy+mv) r(% + mv),i«) ] ^fok)P(y) dy 

, ' n meZ 2 

The above sum can be re-written via the Poisson summation formula as 
^ e -*( k - K My+ mv ) r(£(y + mv), 6s) 



me I 



m\k\ //'_>kj t- ! k. K) 
This completes the proof of Proposition 17.1 



g-2 £ e im. y f [ I^l^g^±^ >Sa ) (7.11) 

m— (mi ,m2)e^ 2 



7.1. Estimation of t)||i, 2 (R 2 ). 

In this section we prove that for any fixed p > and ex > 0, 



sup ||//,z 5 (-^)||l 2 ( R 2 ) = 0{5~), no (7.12) 



Since 



/i,2,(x,t) = 2 / X (|k-K|<^) f b (k,t) $ 6 (x;k) dk, 



we have 



ll//,i?M)lli 2 ( R 2 ) = S / x(|k-K| < sA\f b (k,t)\ 2 dk (7.13) 



be{ + .-} 

where, by (|7.3j) 



/+(k,i) = -z5 2 / ds e -i(n±M-»*)(t-s) J]<$ ± (.;k),7 r («.,«a)* r (.)> 
Jo r=1 



! ) 



(7.14) 

We next use Proposition 17. II to re-express the inner products appearing in (|7.14l) 
as follows: 

<*±(-;k)>7r(M*)*r(0>£a(Ra) 

= jf [ £ 7r ( mikl + W2k / + (k - K) ,^) ] i^tak) P r (y) dy 

(7.15) 
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Since our goal here is to estimate fi ±, we recall that k is restricted to:{|k— K| sS <5 r }. 
Thus, we rewrite the sum f|7. 15[) in terms of its m = and contributions: 

2 (H 2 ) 



1* ^( k - K ) 

Jr 



S 2 



■Ssj P+(y;k)P r (y) dy + j ^E T>5 (|,*«;k) p+(y;k) V r {y) dy 



Termi r 



(7.16) 



where 



<T 2 £ rj (rV,fc;k) 

meZ 2 \{(0,0)} 



4m . y . , miki + m 2 k 2 + (k - K) 
e [ r , OS 



(7.17) 



We now study the terms Terni2 jr and Ttermi r in (|7.16p . 

Estimation of Terni2 ir of (|7. 16[) : We begin with the following 

Proposition 7.2. Denote by j r any of the functions dxctj, dx t C(j, j,l = 1,2. 
For A > 0, assume cto(X) = a(X, 0) e W^ +1,1 (R 2 ). Then, there exist positive 
constants C\,C% such that the following holds: 

1. For any A > and all eel 2 : 

C 

\%(k,T)\ < - r-rp ||a |vF^+i>i(R 2 )' ( 7 - 18 ) 

2. For all A > 2 and a/Z he Bh such that |k — K| < 5 T we have 

| £ r , 5 (y,T;k) | < C 2 <5 A ||a |^ +1 , I(R2) (7.19) 

Proof of Proposition \ 7. 2\ Each function j r itself is the component of a solution 
of the Dirac system. Hence, by Proposition 14. II and integration by parts: 



C 



|7r(«,T)| < C (d x a){K,T) 

J e iK X <5 x a (X) dX 



(dxao) (n) 
1 



< 



(i + M 



-j \\ao\\w A + 1 . 1 (m 2 ) (7.20) 



This proves ([7T8|) . 

We now turn our attention to the bound (|7. 19|) . Note that the sum in the defini- 
tion of E r j, (|7.17[) . is over m = (mi, 7712) £ Z 2 \{(0, 0)}. For such m, and for k such 
that |k — K| < S T , we have 

|miki + m 2 k 2 + (k — K)| 5= c|m| . 

Thus, by Proposition 17.21 

keB ft , k - K < <T => 



I E r , s (y,T;k) I < 



tieZ 2 \{(0,0)} 



wn-y -> , miki + m 2 k 2 + (k- K) 



(7.21) 
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where we take A > 2 for the sum to converge. This completes the proof of the bound 
([TH?)) and therewith Proposition [7^1 

By Proposition 17.21 and (j6.24j) we obtain the following bound on the last term in 
(17161) ': 



k6B ft , |k-K|<<T 



|Term 2 J 



/ S- 2 E r . s (.TV, Ss; k) p+(y; k) V r (y) dy 
Jn 



in 

cA-2 



< 5 A ~ Z \\a \\ wA+1A(W) , A> 2 . (7.22) 



Expansion and estimation of Ternii r of (|T. 16[) : Termi, r may be written as: 

('))l 2 (^) • 



Tferm 1>r = r 2 X (|K-k|<^)-7r( ) • <P±(-;k),P r (-)>2 



(7.23) 



Note that the argument of % may be small and hence we may no longer use the 
decay properties of j r to control the magnitude of (|7.23l) . We therefore use the precise 
behavior of p+(x;k) for k near K. By (P3) 

, ^ 1 (ki-KQ + ifo-Ka) 1 

p±( x ; k = -7^77; — ttt; — — „ , 9 ,i pi( x ± -^P2(x) 

V2 ((ki - Ki) 2 + (k 2 - K 2 ) 2 ) 2 V2 
+ Ojppp/Ah) ( ((k x - KO 2 + (k a - K 2 ) 2 ) ' ) (7.24) 

Next we use (j7.24jl to expand the inner product in (|7.23[) . We have, recalling that 
k = K + k: 



^ ( k 8 K ' Ss ) ' ^±(' ;k )' Pr (')>-t 2 (n) 



1 Ki+iKl / ^ f K \ ^ \ / p ^ ( tl Ss) ■ V r / 

+ C(V ' ll a o||w 2 > 1 (E 2 )) 

1 '■•'-%/<^(" fa ).*A. ±i/*,,5,f 



+ • ||ao||w 2 ,i(M 2 )) (7.25) 

The inner products in (|7.25l) can be evaluated using the expressions for 7 r (X, T), 1 
r ^ 4 displayed in (j6. 16|) - (16. 19[) and the following mild generalization of Proposition 
4.1 of [10] to the case of complex ( = (Ci, C2): 
Proposition 7.3. 

<$a,C- V$ Q >= 0, a =1,2, 



2i <$i,C-V$2>= 2 4 <$ 2 ,C-V$!> = -A» x (Ci+tCa), 
2z<$ 2 ,C-V$ 1 > = -Aj x (Ci-iCa) (7-26) 
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where 

Aj = 3 area(ft) x £ c(m,^) 2 ( 1 ) • K m . (7.27) 



me 5 



Here, {c(m; /i*)}meS denotes the sequence of Fourier coefficients of the normal- 
ized eigenstate 

*i 6 iic.j n Nullspacef - A + Vh - A**J 

and cScZ 2 is as in 

The inner products in (|7.25|) can be evaluated, thanks to (|6.16j) - (|6.19j) and Proposition 

Inner products of the form: <($i , •) : 

<$i,fi • *i> L 2 (n) = id T ai, 
(*i.TS -*2> L 2 (n) = 0, 



<*i,Ta-*3>i»(n) = ($i,2V x ai- V x $i^ = 0, 



<$i,74 • *4> L 2 (n) = ^i,2V x a 2 'V x $ 2 ^ 2(si) = iA B (d Xl a 2 + i d X2 a 2 ) 



(7.28) 



Inner products of the form: ($2 , •) : 

<$2,7l'*l>L2 (n) = 

<$2,72 ' *2> L 2 (n) = i5 T a2, 

<$ 2 ,73 • *3> L 2 (n) = ^$2,2Vxai • V x $x/ = ( 8 *i a i - * <3 X2 ai 

<$ 2 ,74-* 4 > L 2 (n) = (* 2 ,2V^ai-Vx$2) = 0. (7.29) 

In (|7.28p and (|7.29[) the Fourier transforms are evaluated at (j, Ss). 
Now summing both sides of (|7.25p over 1 r < 4 we obtain: 

4 

2 M 7'^) ' <p+( , ' k )' :P '-(')>i 2 (o) 



i «i + «ki 



d T ai + Aj ^Jdtt2 + i dx 2 a2^ ](^ (5s 



(/s? + «i)3 

i K\ + in\ ^ /k \ i /k . \ /!-> / 1 A 



(7.30) 
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Remark 7.1. Note that the Dirac equations (|1.5[) - (|1.6[) are equivalent to the 
equations 

Pi(X,T) = 0, X> 2 (X,T) = . (7.31) 
Recalling the definition (|7.16[) of Ternii. r , we find that: 

4 4 

2 Ter mi , r = • k),P r (-)W) X (|k-K|<<T 

"'-^■^^^[^(^ ■«(?«■) 

+ 0(5- 2 \K\x<\K\<5 r )'Mw>*<*>)) ■ (7-32) 
By (JESJ, ([71i]l . (|7IBjl . (gl2| and (021 . for any fixed i > 0: 

io 

x ^ <*±(-;k),7r(«-,fa)*r(-)>ia(R>) x(|k-K|<« T ) 







-»(M+(k)-M*)(i-s) 

+ o(|t| |k-K| X (|k-K| <5 T )-||ao|| w2 ,i (K 2 ) ) 

+ o(]t] X (|k-K|<<T) i A • OoollwrA+i.i ) (7.33) 
Recall that 



ll//,fl(t)||l 2(R2) = 2 / x(|k-K| <SA \f b (k,t)\ 2 dk 

b=± 

= V f \fj, b (k,t)\ 2 dk (7.34) 

6=+ J ^H 

Remark 7.2. In this remark we assess the contribution of the first term on the 
right hand side of the upper bound (|7.33[) to || /j,£>(t) || x,2^]R2) for times t as p (5~ 1-eo . 
and argue that the Dirac equations (jl.5p - (|1.6p are necessary for f s (x,t) to be small 
on a time scale of the order S^ 1 . 

This contribution to the right hand side of (|7.34j) is bounded by a constant multiple 

of 

\t\ 2 J x(|K-k| <( r) dk * ( P r 1 - 60 ) 2 s 2t * (s-^+^y 

Recall that < r < 1 and therefore this contribution diverges as S J, 0. Indeed this is 
the case even for t ^ S^ 1 , (eo = 0). We conclude that for f to be small in L 2 (M. 2 ) for 
5 I we must have 

x {\k\<S t - 1 ) [tfrai + (d^2 + ifa~a 2 ^ ] (|> 5 *) = 

xQ/e^tf 1 - 1 ) + A tt (^"51 - * fo~ai) ](%> 5t ) = °- ( 7 - 35 ) 
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Indeed, these are implied by the Dirac equations (|7.31l) or equivalently (|1.5[) - (ll.6[) . 

Since the Dirac equations are assumed to hold, /j,±(k, t) is controlled by the 
latter two terms in (|7.33|) . Their contributions to ||//,d(£)||^2( R 2) are bounded, for 
^ t ^ p 6~ 2+£l , as follows. Fix E\ > 0. The second term in (|7.33p gives a 
contribution to ||//,£>(£)||£2(r2) which is bounded by: 

\t\ 2 f |k-K| 2 x (|k-K|<<r) dk 

* \t\ 2 S iT ^ pei-M.i-r) <g S ei ( 7 36 ) 

by taking r chosen sufficiently close to 1. The third term in (|7.33p . for A sufficiently 
large, clearly gives a contribution to ||//,D(i)li2( R 2) which is 0(S £l ) for t 
p 5- 2+£ \ ei > as 6 i 0. 

In conclusion, for any fixed p > and £i > , we have 

sup ||/r,fl(t)||i»(B») = 0(<S^) ■ (7.37) 

This proves the bound (|7.12l) . 

7.2. Estimation of |//z,z)(t)|i,2(] R 2). 

In this section we prove that for any fixed p > and E\ > 

sup ||///,c(-,*)|U3(R3) = as 5i0 (7.38) 

0sStsSp5- 2 + e i 

Since 

/j/, D (x,t) = J] f x(|k-K| > 8 T ) f b (k,t) $ 6 (x;k) dk , 

we have 

ll/i/,D(-,*)lW a ) = E / x(|k-K|>^) |A(M)| 2 dk (7.39) 
where we recall again, by (|7.3p . that 

r* 4 

/±(M) = -*<5 2 / dse-ifeW-^JM ^ <3 ± (.. k ), 7r (5-,5 s )* r (-)> i2(E2) • 
Jo r=1 

(7.40) 

By Proposition 17. II with the choices: 6 = + or 6 = — ,T = 7 r (<5x, (5s) and V = V r we 
have 

<$+(•; k), 7r (5-,5s)* r (-)> L2(R2) 

= jT ^k) [ r 2 £ e— % ^A + TO2 k 2 + ( k-K) ^ ^ Pr(y) rfy 

(7.41) 
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The next lemma implies that all terms in the infinite sum within (|7.41[) involve % 
evaluated at a large argument. 

Lemma 7.4. Let he Bh and assume |k — K| 5j 6 T . Then, there exists a constant, 
c > 0, such that for all m = (mi, mi) e Z 2 

| mxki + m 2 k 2 + (k - K) | > c<5 T (l+ |m|) . 

Proof of Lemma \7.4\ \mik.\ + m 2 ^-2 + (k — K) | is equal to the distance from the point 
K — m]ki — m 2 \i.2 to k. Simple geometry concludes the proof. 

By Proposition 17.21 and Lemma T7. 41 . for k e Bh and |k — K| 5= 8 T we have 



S ~ 2 h < 

meZ 



L. y « / miki + TO2k2 + (k — K) 



<: 



c S n + i m hA ^ 2 ( 5l ~ T ) A II"o||wa+i,i(r2) ■ 

From ([7^40]) and (fT3T]l . we have that for k e B h and |k - K| 5= <T by taking A > 2: 

/±(k,t)| < C |t| laolwA+i.!^) • (7.42) 

By ([7351) . for * < <5~ 2+£l , by taking A > A > 2 sufficiently large, we have the 
bound !///,+ (•, t)lz, 2 (R 2 ) = 0(S 10 ), from which (jT^g]) follows. Recalling (|7T57|) . we 
obtain 



sup ||/r»(t)||i=(M=) = 0(5 2 ) 
0sgts;i5- 2 + E i 



(7.43) 



7.3. Estimation of ||/d<=(-, i)||jja( E 2). 

Again, we fix p > and £\ > and assume < t pd~ 2+£l . Recall from (|7.5 



that 



f D e(x,t) = 2 / A(M) *b(x;k) dk , 
**{+,-} Sh 

/ b (k,t) = ^ 2 



(7.44) 



(7.45) 



To prove that 



SUP ||/d-(-, *)||j2-=(]R2) 

we shall decompose /cc into its k- components near and away from K. For k near 
K, Property (P2), (|6.2I) . implies that the complex exponential in (|7.45[) is oscillatory 
and so integration by parts gains us smallness via additional powers of S. For k e Bh 
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but away from K we use that the Fourier transform of 7 r (<5x, <5s)^E' r (x) is mainly 
supported away from K. Finally, smoothness of 7 r (<5x, <Js)\Ev(x) is used to ensure 
sufficient decay as b — ► oo, of its $&(•; k)- components, which must be summed in 
order to control H s norms. 

To implement the above strategy we decompose fn<= as 



/ D c(x,t) = J] f X (|k-K| < k x ) f b (k,t) $ b (x;k) dk 

m+-} Bh 

+ E / x(|k-K|>«i) A(k,t) $ fe (x;k) dk 
¥{+-} Bh 

= /j,D°(x,t) + //I,D<=(x,t) 

By property (P2), (|6.2p . we have 

|k-K|<«i |/i 6 (k) - /x*| > Ci for 6 ^ ±. 



(7.46) 



(7.47) 



It is natural to exploit the oscillation coming from the complex exponential in (I7.45[) . 
Fix keBj, and such that |k — K| < k%. By (|7.47[) , we may integrate by parts once 
and obtain: 



MM) = s 



■z e 



-»(f*&(10— /**)* 



/Ui,(k) — jU* 



<M-;k),£ 7r («5,0) *r(- 



+ i5* ds 
Jo 



>(Mb( k )-M*)(*-s) _ ^ 



L 2 (R 2 ) 



/i6(k) - /x, 



4 6 (.;k),2 3r>(<M«) * P (. 



t 2 ( 



(7.48) 



Therefore 



|/i(M)| < <5 2 



iMfe(k) -M*| 



*i(.;k), j 7r (*- > 0) * r (0 



1 



+ S 3 I ds 

/o lMfe(k)-M*l 



r=l / L 2 (B 2 ) 

/ 4 \ 

$ 6 (-;k),£dr7r(M«) * r (.) 



L 2 (R 2 ) 

(7.49) 



For |k - K| < ki and b ^ +, we have by (|7.47|) 



l/ b (k,t)| <5 2 



$ b (-;k),27r(^,0) * r (.) 



+ <5 3 |t| max 

0^s^pi5- 2 + e o 



= 1 / L 2 (R 2 ) 

/ 4 \ 

$ b (-;k),2^T7r(^,^) * r (.) 



(7.50) 



L 2 (R 2 ) 
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Therefore, for |k — K| < Ki and b ^ + 



sup |/ 6 (k,t)| 2 



< S 4 



* 6 (.;k),27r(«-,0) %■(■) 



+ 5 2+2eo max 



L 2 (R 2 ) 
4 



$ 6 (-;k),5]Sr7r(^^) *r(0 



r=l 



L 2 (R 2 ) 



(7.51) 



To obtain a bound on ||/7,d c (', ^)||_fr s (R 2 ) for any s > 0, we proceed as follows. 
Recall 

WfiM;t)\\H* m ~ £ (l+H) S /r . IA(k,i)| 2 rfk (7.52) 



^{ + ,-} 

By £/ie sum over 6 ^ { + , — } we mean the sum over all b ^ 1 swc/i i/iai 6 ^ {&i, b\ + 1}, 
where /^(k) = A*-(k) and /i(, 1 + i(k) = /i + (k); see Definition \S.l\ 

Thus we'll require decay of |/&(k, t) \ for & large. This decay is obtained from the inner 
products in (|7.51l) . Observe, for some sufficiently large postive constant, C and any 
M 5s and j S* 0: 

*6(-;k),4>(*-,**) *r(0 
1 



i 



(C + Mf ,(k)) M 
Therefore, thanks to (pT2Tj) : 

$ 6 (-;k),^ T 7 r (5-,5s) * r (.) 



L 2 (R 2 ) 

(C + J ff) M <E- b (-;k),^7r(5-,5 S ) * P (.) 
$ 6 (-;k),(C* + //) M f4>(«.,tfa) * r (.) 



L 2 (R 2 ) 



L 2 (M 2 ) 



(l + |b|) _M tr(8;5s) * r (.) || 



< (1+|&[) M S 1 ||ao||H2M + j( R 2) 



for j Ss 0. Using ([733)1 in ([73T1) we obtain 

sup \f b (k,t)\ 2 < (l+|fe|)- 2M \\a \\ 2 H2M+1(R2) S 2 ^ 

OsStsSp<5- 2 + E o 

Substituting into (|7.52|l . we get 

WiM-Mb-OP) ~ ( S (l + |6|r 2M )' ||aob- 

&*{+-} 

< ||ao||_H-2Af + i( R 2-, 5 £ ° 

provided 2M > s + 1. 

It remains to estimate the H s norm of 

x(\k-K\^Kx) A(k,t) $ 6 (x;k) dk, 



(7.53) 
(7.54) 



(7.55) 



+ ,-} Bh 



(7.56) 
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where we recall 

\\fu,Dc(-;t)\\ H , 



2 (l + \b\Y [ x (|k-K|> Kl ) \f b (k,t)\ 2 dk . (7.57) 



b<f{ + .-} 



Note that from §TM\i 

X ( |k-K| =s «i) |/ 6 (M)j 2 dk 



s£ <T|tr max 

Osjssjt 



X |k-K| > m 



L 2 (R 2 ) 



dk (7.58) 



A crude bound on (|7.58[) . valid for < t < <5- 2 + £ o uses the approach taken to 
obtain (|7.53|) . This gives the bound: S 4 (5~ 2+e ") 2 5~ 2 ||ao|^a,i(Ha)) which becomes 
unbounded as (5 J, 0. Therefore, a sharper estimate is required. 

Now, for k e B, |K— k| > K\ and 6 ^ +, it may be that |/i(,(k)— //*| is small. Hence 
the integral in (I7.45|) cannot be controlled as an oscillatory integral, via integration 
by parts with respect to time. We therefore obtain the decay of fii,D c using the rapid 
decay of the Fourier transform of "f r . 

The inner product can be rewritten and estimated as follows. Choose a positive 
constant C such that CI + H is strictly positive. We have, for any integer M > 0, 

<$;,(•; k), 7r (<5-, 8s) ■ 4v(-)>l 2 (r 2 ) 



(C + ^(k))M 
1 



_L 2 (R 2 ) 



* 6 (-;k),(C + H) A1 [ lr (8;6s) * r (.) 



\L 2 (R 2 ) 



(C + ^(k))M 

where J] denotes a finite sum over terms of the above form with a\ , . . . , a„ max , b\ , &2 
in and 2^ maa; + 2^™^ |a„| + + I&2I 2M. Each inner product of this sum can 
be re-expressed, via Poisson summation, using Proposition 17.11 With the choices 

r(X,T) =^ 7r (£x,<5 S ), 

*(x) = n^r^^(x) S^* r (x) = e iKx n^fWx) + *K) b ^,,(x) e L 2 K 
we apply Proposition 17. II and obtain 



/ Pft (y;k). [<r 2 ^ e toy d£ 



7r 



L 2 (R 2 ) 

miki + ra2k2 + (k — K) 



(7.60) 



Note that there is a constant c > 0, depending only on k\, such that for all k e Bt 
satisfying |k — K| 5= «i, we have |miki + m2k2 + (k — K)| 5= c(l + |m|). It follows 
from Proposition 17. 21 applied to (I7.60[) that for all b\ and all A such that A > M and 
A>2: 



- 5/1 ll"o||w-^+i.i(R 2 ) I] (1 + |ml1 A ~ ^ H Q °II H/A+1 ' 1 ( 

meZ 2 ^ ' " 



L 2 (R 2 ) 



(7.61) 
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Recall that /ifc(k) sb b, for b large, uniformly in k e £?/,. Using (|7.59|) and (|7.61l) 
in (f737j) and (1735)) we obtain 



-2M 

V-l -r 



1/77,^(^)1^ S ^ Kl 2 ^ l«0|^ +1 , 1(R2) 2(l + &) 

Choosing A and M sufficiently large, we obtain 

sup \\f IItD c(-;t)\\ Hs < S 10 \\a \\ W A + i,i(R2) . (7.62) 

Together, estimates (|7.55[) and (|7.62[) imply 

sup \\fD4-,t)\\Hs ^ 5 £o \\a \\ W A+i,i( R 2) . (7.63) 

Finally, ([7]35]l . (17331 and (j2~2"5|) imply, for any e > 0, p > 0, that 

sup ll/' 5 (*)l|77«(B2) ^ ^ ||ao||w^+ 1 - 1 (K 2 ) ( 7 - 64 ) 

OsS*sSp(5- 2+e o 

This completes the proof of Proposition 16.11 and therewith Theorem 15.11 

Appendix A. Lipschitz-continuity of eigenvalues self-adjoint 2 nd order 
elliptic operators and an application to Floquet-Bloch eigenvalues. 

Theorem A.l. Let T and T denote operators with the following properies: 

1. T and T non-negative and self-adjoint operators on L 2 . 

2. T and T are bounded maps from H 2 to L 2 . 

3. T and T have discrete spectrum given, respectively, by the sequences of eigen- 
values: 

spec(T) : A X (T) < A 2 (T) < • • • 
spec(T) : X 1 (T ) < A 2 (T) < • • • 

(A.1) 

4- There is a positive constant, , such that T satisfies the elliptic estimate: 

\\Tv\\ L2 + \\ip\\ L 2 >a%\\ H i (A.2) 
for all (p e H 2 . Assume that 

\\T-f\\ H 2^ L 2^lj (A.3) 
Then, for k = 1, 2, . . . we have the following Lipschitz estimate 

X k (T)-X k (f) *S -{ 2A fe (T) + l )-\\f-T\\ H ^ L 2 (A.4) 
a* 

Remark A.l. Theorem I A. 1\ generalizes in a straightforward manner to higher 
order self-adjoint elliptic operators, defined on H 3 (M. d ). 

A consequence of Theorem IA.1I is the Lipschitz continuity of the Floquet-Bloch 
eigenvalues of H = —A + V(x), where V is periodic, real and bounded. 



32 Effective Dirac Dynamics 

Corollary A. 2 ( Proposition 12. IT) . The Floquet-Bloch eigenvalue maps k — > 
/i{,(k), 6^1 are Lipschitz continuous functions ofkeB. 

Proof of Corollary \A.2c Define 

T = Jf(ki) = -(V + iki) 2 + 7(x) and T = ff(k 2 ) = -(V + ik 2 ) 2 + V(x) 

and note that \\T-f\\ H *^, L 2 = \\H(ki) -H(k 2 )\\ H ^L^ < C 1 |ki-k 2 | for some C > 0. 
Let /x&(kj) = fib{H(kj)) denote the eigenvalue of H(kj). Applying Theorem lA.il 
we have 

Mkx) - w (k 2 )| *S C ( \» b (k)\ + 1 ) |k x - k a | (A.5) 

whenever |ki — k 2 | is less than a small enough positive constant. This completes the 
proof of Corollarv lA.2l 

Proof of Theorem \A.1[ Note first that 

II^VIl 2 + IMU 2 > [ ||^|| L 2 + \\lfi\\ L 2 ] - \\T -f\\ H 2^ L 2\\(p\\ H 2 

>a%\ B? - \\T-f\\ H ^ L 4<p\\ H 2^^\\y\\ H 2 . (A.6) 
Therefore, we have the following three estimates: 

IIVI^ + MU» > \aHv\\ H * (A.7) 

||2V|U« + \\ V \\ L 2 > ^a?\\<p\\ H 2 (A.8) 

\\T-f\\ H 2^ L 2^lj (A.9) 

These conditions, which are symmetric in T and T, will be used below. 

Recall now the min-max characterization of the k th eigenvalue, Xk(A), of a self- 
adjoint operator, A : H 2 — * L 2 [5]: 

AfcW) = min max (A. 10) 

Scff 2 , dim(S) = k t>eS\{0} ||-u||^ 2 

Proposition A. 3. Let j4 denote a self-adjoint operator which maps H 2 to L 2 
satisfying the ellipticity estimate (|A.7|) . 

X k (A) = min max ( A - n ) 

Scff 2 , djm(S) = fe ¥>eS\{0} ||vIl2 

for all -ipeS 

Proof of ' Provosition \A~3c First note that the min max in ([A. Ill) , greater than or equal 
to Afc(.A), given by (|A.10[) . We claim A& is achieved at an eigenfunction, ipk ¥= 0, with 

ijjk 6 Sk(A) = span of the first k cigenfunctions of A, 

where any tp e Sk(A) satisfies: 

\m H 2^^)-\x k (A) + im L 2 . (A.i2) 



Charles L. Fefferman and Michael I. Weinstein 



33 



Indeed, for any tp e S k (A), the span of the first fc— eigenfunctions of A, we have 

\\Aip\\ L z < A fe (A) \m L2 . 

It follows from (|A.7j) that 
3 



£ 2 



< (A fe (A) + 1) IVII 



L 2 • 



Thus, any ?/> e Sk(A) satisfies (|A.12[) . Furthermore, the maximum of the quotient 
{Aip, ijj) /I V'lia over <5'fc(^)\{0} is equal to Xk(A) and is attained at ipk- This completes 
the proof of Proposition IA. 31 

Continuing with the proof of Theorem lA.il take 5* to be any subspace of dimension 
fc and such that 



veS = 
For all ^ v e S, we have 



\\v\\ m < -(aT^Afc^ + lJIMU,. 



Therefore, 



(Tv,v) L 2 




V w 'Vl 2 


Ml 2 














Ml, 




:> 





















+ 



V- 



L 2 



11^- T|| H 2^ L 2 ||tj|| ff 2 ||tl|| i2 



A fc (r) > A fe (f) 



L 2 



-|r-r|^ £ . ^- 



- \\T - f\\ H ^ L2 ■ -(a*)- x (\ k (T) + 1) (A.13) 



|T - f || ff2 ^ L2 • -(o'j-^AfcCT) + 1) , (A.14) 



thanks to (jA.llj) with A = T and ([A.10I) with A = f. Interchanging f and T in the 
above argument yields 



A fc (T) 5* A fe (T) -iT-rlfla^.-CaT^AfcCTj + l) 



(A.15) 



Estimates (|A.14j) and (|A.15|) imply that 



A fe (T) - A fe (f) *S -{a»)-\\ k (T) + X k {f) + l)-\\T-f\\ H ^ L 2 



< -(a»)^ 1 (2A fc (T) + l).|T-f||^^ 2 



+ -(a")- 1 -\\T-f \\ H ^ L2 



A fc (T) - A fe (T) (A.16) 



Therefore, since we have assumed IIT — T\\h 2 ^l 2 ^ la , we find that 



A fe (T) - A fe (f) < -(a»)- 1 (2A fc (T) + 1) • ||T - f \\ H ^ L2 
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This completes the proof of the Lipschitz bound (|A.4[) and therewith Proposition lA.il 
Finally, Proposition 12.11 is an immediate consequence of Proposition IA.1I 
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